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ABSTRACT 

We investigate the short-term dynamical evolution of stellar grand-design spiral arms 
in barred spiral galaxies using a three-dimensional (3D) A-body/hydrodynamic sim¬ 
ulation. Similar to previous numerical simulations of unbarred, multiple-arm spirals, 
we find that grand-design spiral arms in barred galaxies are not stationary, but rather 
dynamic. This means that the amplitudes, pitch angles, and rotational frequencies of 
the spiral arms are not constant, but change within a few hundred million years (i.e. 
the typical rotational period of a galaxy). We also find that the clear grand-design spi¬ 
rals in barred galaxies appear only when the spirals connect with the ends of the bar. 
Furthermore, we find that the short-term behaviour of spiral arms in the outer regions 
(i? > 1.5-2 bar radius) can be explained by the swing amplification theory and that 
the effects of the bar are not negligible in the inner regions (i? < 1.5-2 bar radius). 
These results suggest that, although grand-design spiral arms in barred galaxies are 
affected by the stellar bar, the grand-design spiral arms essentially originate not as 
bar-driven stationary density waves, but rather as self-excited dynamic patterns. We 
imply that a rigidly rotating grand-design spiral could not be a reasonable dynamical 
model for investigating gas flows and cloud formation even in barred spiral galaxies. 

Key words: method: numerical — galaxies: kinematics and dynamics — galaxies: 
spiral — galaxies: structure 


1 INTRODUCTION 


Spiral arms in disc galaxies can be classified into three cat- 
egories: grand-design, multiple -arm, and flocculent spirals 
llElmegreen fc Elmeereenl Il987l ). Grand-design spiral arms 
are large-scale coherent, symmetric, two-armed patterns 
(e.g. NGC 628), whose number fraction is more than 50% 
in ne arby spiral galaxies llGrosbol et al.l[2004l : iKendall et al.l 
Observations show that gran d-design spirals are as¬ 
sociated with bars or companions llK ormend v &: NormanI 
I 1979 I : ISeigar fc Jamedfiggsl : [Kendall et al.ll201ll ). suggesting 
that a bar or companion is essential in forming a grand- 
design spiral in a disc galaxy. In barred spiral galaxies, 
the grand-design spirals extend from the ends of the bars. 
Furthermore, although some studies found little or no cor- 
relation between bar strengths and spira l arm strengths 
dPurbala et al.l l2009l : iKendall et al.| 201 ll). other s tudies 
have found correlations ( Block et al.|[2004l : [salo et al.ll201^ . 


These observations, therefore, suggest that bars drive grand- 
design spirals with the same pattern speed as the bars. 

This ‘bar-driven spiral hypothesis’ has been promoted 
by the ballistic closed-orbit theory, which is based on non- 
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self-gravitating hydrodynamic simulations in fixed barred 


potentials (e.g. Sanders fc Huntl^ 19761 : ISandersI Il977l : 


iHuntlev et al.l1l978l : Sormani et al.l l2015l ). This theory ex¬ 


plains the physical mechanism responsible for the formation 
of the spirals in b arred galaxie s in terms of closed orbits 
of the gas (see also I Wada 1 19941 ) : the spiral arms in barred 
galaxies are regarded as kinematic density waves (i.e. crowd¬ 
ing of gaseous closed orbits) driven by an external barred 
potential. However, spiral arms in observed barred galaxies 
are composed of not gas but stars; therefore, theories based 
on stellar dynamics are required for investigating the origin 
of spiral arms in barred galaxies. 


On the other hand, the ‘inv ariant manifold theory’ is an¬ 
other bar-driven spiral theory | Romero -G6mez et al.l l2006l . 

M 


l2007l : [Athanassoula et al.ll2009al i^. I 2 OI 0 I ) and has been devel- 

oped from ste llar orbital theories in fixed barred potentials 
dDanbvIl 19651) . Essentially, the back-bones of barred spirals 
are bunches of untrapped stars escaped from unstable La- 
grangian points Li and L 2 close to the ends of the bar. 
This means that stars should move along the arms rather 
than across the arms (lAthanassoulal I 2 OI 2 I ). This behaviour 
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is completely different from the quasi-stationarj0 density 
wave theory, which predicts that stars should traverse the 
arms b ut would stay lon ger In the arm than in the inter-arm 
region dLin fc Shulll964l i. Furthermore, the manifold theory 
predicts a rigidly rotating spiral arm with a similar pat¬ 
tern speed to that of the bar (i.e. flgpirai — f^bar) and that 
stronger bars should have more open sp irals comp ared to 
weake r bars llAthanassoula et alj l2009bf l. Recently, IStruckI 
(l2015ll extended the classical epicycli c orbit appro ximation 
to the ‘p-ellipse’ orbit approximation (IStruckl2006l ) and pro¬ 
posed that ensembles of eccentric resonant orbits excited in 
Lindblad zones can provide a backbone for g enerating kine - 
matic bars and spiral waves. Furthermore, IStruc 5 12 OI 5 II 
showed that some resonant eccentric orbits have a similar 
appearance to invariant manifold orbits. 

The manifold theory and eccentric resonant orbits the¬ 
ory, however, do not consider the self-gravity of stars, despite 
the fact that the self-gravity of stellar discs plays an essential 
role in the dynamical evolution o f spirals via the (non-linear 1 
swing amplification m e chani sm (|SellTOO(U^_CMlberg| 1984 
Carlberg_&_^Yeedmm]|_^j985L^hju_etah 201ll : iBaba et ahl 


2013l:ln’Onghia et al.l(2ni3l:lD’OughiJ2ni,4> . The swin ^ ; 


plification mechanism, proposed by Toomre in Il98ll . oper¬ 
ates through a combination of three aspects, i.e. the shear¬ 
ing ffovw epicyclic motions, and the disc self-gravity, in disc 
galaxieqj. Because the am plification is most efficient around 
the corotating radius le.g. 1shr]ll992h . prominent spiral arms 
tend to have differentially rotating patterns, which almost 
follow the galact ic rotation at every radius (flspirai — fl; 


_ LSi _ 

iBaba et ahl I20l4 l. Besides, the pitch angles are expected 
to depend on the shear rates of galactic discs (lFuchsll200ll : 
iBaba et al.ll20l4 iMichikoshi fc Kokubol[20l4l . 

For investigating the origin of grand-design spirals 
in barred galaxies, A^-body and A'-body/hydrodynamic 
simulations are required, and several preliminary numer¬ 
ical studies have been proposed. However, the studies 
have not reached a consensus about the dynamics of 
barred spirals. Some numerical simulations of barred spi¬ 
ral galaxies have shown that the bar and spiral arm are 
rigidly rotating patterns with dif ferent pattern speeds, and 
they are independent p atterns (ISellwood fc Snarkel 1 19881 : 
iRautiainen fc S alolll99 9ll or are coupled via non-linear in- 
teractions i Masset fc Taggeilll997l : IRautiainen fc Saldll999l : 


^ In general, the word stationary or steady means ‘not moving.’ 
However, in this context, this word indicates that density waves 
do not propagate radially and that instea d they propag a te az - 
imuthally with a single pattern speed. See iBertin &: LinI lll996l l 
for further discussion on the concept of quasi-stationary density 
waves. 

^ The term ‘swing amplification’ is used to express the ampli¬ 
fication process, which combines the shearing flow, epicyclic mo¬ 
tions, and disc self-gravity, regardless of whether the process is 
a linear or nonlinear phenomenon. Originally, the ‘swing amplifi¬ 
cation’ and its feedback cycle were bas ed on the linear perturba¬ 
tion a nalyses of local shearin g sheets jGol dreich &: Lvnden-Belll 
I 1965 I : ljulian fc Toomi4 1 19661 : iToomrel Il981^ .' and previous N- 
body s imulations of discs partly supported this mech¬ 

anism llCarlberg fc Freedmanl [TqSSI : iFuiii et ahl [201 ih . However, 
A^-body simulations of galactic discs suggested the importance of 
nonli nearity, such as mode-mode coupling and radial migration of 
stars teaba et al. I I 2 OI 3 I: iD’Onghia etaDr2013r> . in the dynamical 
evolution of spiral arms. 


iMinchev et al.l l2012l'l . Recently, iRoca-Fabreea et al.l (l2013fl 
performed V-body simulations of pure stellar barred spi¬ 
ral galaxies and analysed the spiral rotation frequency 
from selected snapshots where the spiral amplitude is 
above 70% of the maximum amplitude; they concluded 
that the spiral rotation frequency approaches the bar’s 
rigid body rotation in barred galaxies with the incre¬ 
ment of bar strength. On the other hand, iBaba et al.1 
(l2009ll performed V-body/hydrodynamic simulations of a 
barred spiral galaxy and suggested that the spiral arms 
are not rigidly rotating patterns but transient recurrent 
patterns (i.e. dynamic spirals), and these results are sim- 
ilar to those obtained from multiple arm spiral galaxies 


llSellwood & Garlberd Il984 iFuiii et al.1 201 ll: Wada et al. 

2011 

: iGrand et al. 2012bl: iBaba et al.ll2013l: iD’Onehia et al. 

2013 

1. iGrand et al.1 (|2012al also found that spiral arms in 


barred galaxies are dynamic patterns whose rotation fre¬ 
quencies decrease with the radius in such a way that the 
rotation frequency is similar to the rotation of stars. 

In this study, to address the short-term (e.g. a few 
hundred million years) behaviour of grand-design spirals 
in barred spiral galaxies, we performed a three-dimensional 
(3D) V-body/hydrodynamic simulation of a Milky-Way-like 
barred spiral galaxy. This paper is organized as follows. We 
describe the galaxy model and our numerical simulation 
methodologies in Section [2] In Section [3] we present our 
results concerning the short-term behaviour of grand-design 
spirals in the barred galaxy and effects of the bar on spiral 
dynamics. Finally, we summarize our results in Section [d] 
Long-term (e.g. ~ 10 Gyr) behaviour of spirals In barred 
galaxies, as well as effects of a ‘live’ dark matter (DM) halo, 
will be presented in a forthcoming paper (M. S. Fujii et al. 
in preparation). 


2 NUMERICAL SIMULATIONS AND 
ANALYSIS 

We performed a 3D A-body/hydrodynamic simulation 
of a Milky Way-like galaxy with an A-body/smoothed 
particle hydrodynamics (SPH ) simulation code ASURA-2 
(ISaitoh fc Making l2009l . l2010h . In this study, we focus on 
the dynamical evolution of spiral arms in a simulated barred 
spiral galaxy, while detailed studies on spatial distributions 
of the interstellar medium (ISM) around spiral arms as well 
as the structures of the ISM/clouds and their galactic en¬ 
vironmental dependences are separa tely discussed in Baba, 
Morokuma-Matsu i fc Eg usa ll2015al ) and Baba, Morokuma- 
Matsui & Saitoh (l2015bl L respectively. 


2.1 Numerical methods 

The self-gravities of stars and SPH particles were ca lculated 
by Tr ee with the GRAvity PipE (GRAPE) method (iMakind 
I 1991 II . The Tree method allows us to reduce the number of 
particle pair interactions that must be computed by dividing 
the volume into cubic cells, so that particles in distant cells 
are treated as a single large particle centred at the cell’s cen¬ 
tre of mass. The GRAPE tool is a special-purpose hardware 
for calculating gravitational forces. In this study, we used a 
software emulator of G RAPE, known as Phantom-GRAPE 
llTanikawa et al1l2013l L 
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Radiative cooling and heating, star formation, and stel¬ 
lar feedback from Type II supernovae and Hn regions were 
also included. The radiative cooling of the gas was solved as¬ 
suming an optically thin cooling fu nction for a wide te mper- 
ature range of 20 K < T < 10® K (IWada et al.ll200^') . Star 
formation from cold, dense gas and feedback models was 
based on the simple stellar population (SSP) approximation 
with the Salpeter initial mass fun ction. These were imple- 
mented in a probabilistic manner llSaitoh et al.|[200^ . The 
Hii-region feed back was implemented using a Stromgren vol¬ 
ume approach (iHopkins et al.l[2012l : iBaba et al.ll2015bh . 


2.2 Galaxy model 


We generated the initial axisymmetric model of a stellar 
and gaseous discs with a classical bulge embedded in a DM 
halo. Fig.[T] shows the initial circular velocity curves of each 
component in the model. 

The stellar disc follows an exponential profile: 


where Md, Rd, and Zd are the total mass, scale-length, and 
scale-height of the stellar disc, respectively. Using Hern- 
quist’s method (lHernauistlll993ll . the velocity structure of 
the stellar disc in cylindrical coordinates is approximately 
determined by a Maxwellian approximation. We set the ref¬ 
erence radial velocity dispersion by assuming tha t Toomre’s 
Q at R = 2.5Rd equals 1.3 (see section 2.2.3 of iHernauistI 
Il993h . 

The gas disc also follows an exponential profile with a 
total mass of Md,g, scale-length of Rd,g, and scale-height of 
Zd,g- The gas disc was truncated at 5Rd,g- The initial tem¬ 
perature was set to 10^ K, and the initial velocity dispersion 
was 10 km s“^ everywhere. 

The classical bulge follows the Hernquist profile: 


Pbulge(r) 


Mb,0 CLh 
2tx r{r + at)® ’ 


( 2 ) 


where Mbfi an d at are the tot al mass and scale-length 
of the bulge (iHemciuisf Il99d l. respectively. Following 
IWidrow fc Dubinskil ( 2005 1. we generated the classical bulge 
from the distrib ution function with an en ergy cutoff (see 
Equation (10) of lWidrow fc Dubinski|[2005ll . Here, we set Qb 
of that equation to be 0.21. 

For simplicity, we assumed the DM halo to be a static 
potential, whose density profile follows the Navarro-Frenk- 
White profile: 


Mb 


Phalo(r) + 

fc{x) = ln(l -I- a:) — a:/(l -I- a:). 


(3) 

(4) 


where Mb, ah, and Cu are the total mass, virial radius, 
and conce ntration parameter o f the of the DM halo, re¬ 
spectively dNavarro et al ] fi997h . and X = rjah- As shown 
in Fig. [T] the contribution from the stellar disc to the total 
circular velocity in the galaxy model surpasses that from 
the DM haloes in the regions of 3 < i? < 12 kpc. In 
this case, the galaxy satisfies t he criteria for the bar in¬ 
stability (lEfstathiou et aHll982f) and spontaneously devel¬ 
ops a stellar bar. Note that a static halo omits dynami- 
cal friction on the bar (e.g. IPebattista fc SellwoodI l2000l : 



Figure 1. Initial circular velocity curves of each component of 
the galaxy. 


lAthanassoula fc Misirioti^ l2002fl : however, dynamical fric¬ 
tion will be weak because the central density of the DM 
halo in this study is low. We will investigate effects of a 
‘live’ DM halo on the dynamical evolution of barred spirals 
in a future study. 

The initial numbers of stars and SPH particles were 
6.4 million and 4.5 million, respectively, and the particle 
masses of stars and SPH particles were about 9.1 x 10® Mq 
and 3 x 10® Mq, respectively. We used a gravitational soft¬ 
ening length of 10 pc. Following the above procedures, 
we generated a galaxy model in a near-equilibrium state 
and then subsequently allowed the galaxy to evolve for 
~ 10 — 20 rotational periods under t he constraint of a xisym- 
metry dMcMillan fc DehnenI l2007l : iFuiii et aP I 2 OIIII . This 
equilibrium state was used as the initial condition for the 
numerical simulation. 


2.3 Quantification of spiral arms 

The quantities that characterize the spiral structure are the 
number of spiral arms, their shape (i.e. logarithmic, hy¬ 
perbolic, segments of lines), and their amplitude (or arm- 
to-inter-arm contrast). In order to analyse these quanti¬ 
ties, we used the one-dimensional (ID) Fourier decompo- 
sition method with r espect to the azimuthal direction (see 
iRix fc Zaritskvl[l995h . More specihcally, we divided the stel¬ 
lar surface density /i(R, 0) into strips using a polar coordi¬ 
nates [R, (f)) and then performed a ID Fourier decomposition 
of the each strip with respect to the azimuthal direction: 

BmiR) = j exp[-im4>]dcf), (5) 

where fl{R) is the azimuthally averaged surface stellar den¬ 
sity at radius R, and Bm{R) denotes a complex amplitude 
of the m-th mode. 

Radial profiles of the pitch angle were calculated using 
radial changes of phases of Bm{R)- If we define the phase as 
4>m{R) = Si[ctan{Bm/\Bm\), then we can evaluate the pitch 
angle at radius R using the following equation: 

, . / _ r^^4*'nT-{R) 

C0t2pj7^(R) — R -jrr . 

art 


( 6 ) 
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We analysed the angular speeds of the spiral arms with 
the Fourier density peak method. In this method, the rate 
of angular phase changes of i?m, flphase.mffl, t), is referred 
to the ‘angular phase speed’ ifWada et al] I 2 OIII : iBaba et al.l 
l2013l : [^ca-Fabrega et al.]l2013tl . 


3 RESULTS AND DISCUSSION 

We present the results of the dynamical evolution of barred 
grand-design spirals. Fig. [5] shows the time evolution of the 
stellar mass distribution over a period of 1.0 - 3.4 Gyr. The 
simulated galaxy has an m = 2 grand-design spiral without a 
stellar bar until t ~ 1.2 Gyr, begins to show bar instabilities 
around t ~ 1.3 Gyr, and then reaches a quasi-steady state 
until t ~ 2 Gyr. As shown in previous numerical simulations, 
we confirm that the grand-design spiral is neither steady 
nor always connected with the stellar bar. For example, the 
m — 2 grand-design spiral is clear at t = 2.0 Gyr, becomes 
faint at t = 2.1 Gyr, and appears clearly again at t = 2.2 
Gyr. This unambiguously shows that grand-design spirals 
change their strengths with a period of about 200 Myr. 

In the following subsections, we first describe the short¬ 
term (i.e. < 200 Myr) behaviour of grand-design spirals and 
compare these results with the predictions from the swing 
amplification and invariant manifold theories. 

3.1 Short-term behaviour of grand-design spirals 

The left panels in Fig. [3] shows the time evolution of the 
spiral arms in the polar coordinates {R, (p) from 2.70 - 2.84 
Gyr. Thisspan corresponds to one cycle of a bar-spiral con¬ 
nection event. During this term, the bar and strong grand- 
design spiral appear at t ~ 2.76 - 2.78 Gyr. If we focus on 
the grand-design arm placed at (p = 180° - 270°, we see 
that this arm is initially short and weak {t = 2.7 Gyr) and 
then grows into a grand-design spiral arm with winding until 
t = 2.74 Gyr. 

The right panels in Fig.[3]show the time evolution of ra¬ 
dial profiles of the amplitude (|i3m=2|), pitch angle {im= 2 ), 
and angular phase speed (Dphase,m= 2 ) of the m — 2 mode. 
We observe that, at the clear grand-design spiral phase (i.e. 
t ~ 2.76 - 2.78 Gyr), the spiral parameters are nearly con¬ 
stant with respect to the radius or slowly decrease from 
3 kpc to the outer regions. In contrast, at the damping 
phase (t > 2.8 Gyr), the radial profiles of all spiral pa¬ 
rameters show radial modulations. Note that the radially 
changing pitch angle at the clear barred-spiral phase seems 
to contradict the traditional picture, where the observed 
geometries of spiral arms are frequently fitted by a log¬ 
arithmic spiral (i.e. cons tant pitch angle; iKennicu ttl ll98ll: 
iKennicutt fc Hodgdl 1982 1. However, as per recent measure¬ 
ments of observed spiral galaxies, the pitch angle for most 
gala xies decreases as the distance fro m the centre increases 
(e.g. ISavchenko fc Reshetniko^ l2013l l. which is consistent 
with our results. 

The top row in Fig. [4] shows the time evolution of the 
amplitudes at each radius. As described above, the ampli¬ 
tudes change with a period of about 200 Myr at each radius, 
although the amplitude at i? = 2 kpc (i.e. the bar region) 
is more stable than those of the other regions. The upper 
horizontal axis in each panel measures the time normalized 


by the rotational period Trot at each radius. We see that 
the amplitude at each radius changes in time spans of 1 - 2 
Trot; that is, the grand-design spirals in barred spiral galax¬ 
ies are not quasi-stationary patterns, but transient recurrent 
structures (i.e. dynamic spirals )0- 

As shown in the middle row in Fig. [4] the pitch an¬ 
gles also change similarly, but their phases are inversely 
correlated with the phases of the amplitude changes. In 
other words, spiral arms become strong as they wind from 
open to closed and then become damped after reaching the 
maximum amplitude. Such evolutionary behaviour between 
|Bm= 2 | and im =2 is clearly shown in the bottom panels in 
Fig. [5] which shows the values of |Bm,= 2 | and im =2 at i? = 4, 
6, 8, and 10 kpc every 1 Myr. We see that, at all radii, the 
amplitudes are very small in leading spirals > 90°) 

and tightly winding trailing spirals (im =2 Pp, 10°), while the 
amplitudes are large in open trailing spirals (im =2 — 20° - 
45°). 

The bottom panels in Fig. [4] show the time evolu¬ 
tions of the phase speeds, Hphase,m= 2 , at each radius. The 
phase speeds at i? < 3 kpc are stable and approximately 
45 km s“^ kpc“^ because the bar dominates in these re¬ 
gions. In contrast, the phase speeds both increase and de¬ 
crease compared to the galactic rotation for all radial ranges, 
and as for the evolution of the pitch angles, the phases are 
inversely correlated to the phases of the amplitude changes. 
The bottom panels in Fig. [6] clearly show such evolutionary 
behaviour between the phase speed and amplitude. When 
the phase speed at each radius is approximately equal to 
the galactic angular speed, fl(i?), the amplitude reaches a 
maximum. 

Finally, we discuss the evolution of the bar. As shown 
in the left panels in Fig. [31 the bar with a semi-major axis of 
about 3 kpc is always located around (p = 0 and 180°, even 
though its parameters show the time evolution. The pitch 
angles and phase speeds of the right panel in Fig. [3] show 
flat radial profiles at i? < 2 kpc. In particular, at the clear 
grand-design spiral phase (i.e. t ~ 2.76 - 2.78 Gyr), we can 
see that these flat profiles reach i? ~ 3 kpc. As shown in 
Fig- 111 the amplitude and angular phase speed at i? = 2 kpc 
are relatively stable, suggesting that the inner structure of 
the bar is steady. However, the amplitude and angular phase 
speed clearly change in time at i? = 3 kpc. Such behaviour of 
the stellar bar is usually seen in A-body simulations of iso¬ 
lated discs (iMartinez-Valpuesta et al]l2006l : iDubinski et al.l 
I 2 OO 9 II as well as cosm ological simulations of disc galaxies 
( Okamoto et al.]l2015h . 


3.2 Comparison with swing amplification theory 

We compare our results with the swing amplification the¬ 
ory to investigate the physical origin of the short-term evo¬ 
lutional behaviour among spiral parameters. Before this 
comparison, we briefly review the physical essence of the 


® For unbarred spiral galaxies, ISellwood fc Carlberd ll2014^ ar¬ 
gued that the dynamically evolving spiral patterns are attributed 
to the superposition of a few longer-lived waves, each of which 
has well-defined frequencies and shapes and survives for 5-10 
rotations. 
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Figure 2. Face-on stellar mass distributions of the galaxy every 200 Myr from 1.0 Gyr to 2.0 Gyr, and every 100 Myr from 2.0 Gyr to 
3.4 Gyr. The bar instability occurs around t ~ 1.3 Gyr. After the bar is well developed, grand-design spirals are clear around 2.0, 2.2, 
2.4, 2.6, 2.8, 3.0, 3.2, and 3.4 Gyr. 


swing amplification (see iDobbs fc Babal l2014l . and refer¬ 
ences therein for details). The swing amplihcation operates 
through a combination of the shearing flow, epicyclic mo¬ 
tions, and disc self-gravity, which play roles in promoting 
or s uppressing the formation of structures in disc galax¬ 
ies (iToomrel Il98lli . Since the direction of epicyclic mo¬ 
tion of a star is the same as the direction which the spi¬ 
ral arm is sheared by differential rotation, stabilisation by 
rotation is reduced, and t he perturbation can grow via 
the usual Jeans instab ility llGoldreich fc Lvnden-Belllll965l : 
ljulian fc Toomr^ll966l ~). As a concrete example, we consider 
a disc galaxy with a nearly flat rotation curve as our simu¬ 
lated galaxy (see Fig. [!}. Because the angular speed follows 
as oc R~^, the timescale of epicyclic motion is given by 


_ 1 
T~epi = — 
K 




- 1/2 


1 

v/2n’ 


( 7 ) 


while the timescale of involvement with the spiral arm is 
given by 


R\dn/dR\ n' 

Comparing the above two time-scales, we can see that 
Tepi ~ Tshear- The shearing motion cancels out the stabil¬ 
isation effect of epicyclic motions, and thus, the structure 
can grow via self-gravity in a short time as compared to 

'7”epi- 

To compare the results described in Section 13.11 with 
the predictions from the swing amplification theory, we cal¬ 
culate the pitch angle at the maximum amplitude phase, 
^s wA, expected from the swing a mplification theory. Follow¬ 
ing [Michikoshr&^^oku^ ll2ni4h . iswA is given by 


*swA = arctan 


/2 V4- 2r\ 

V7^^ J 


( 9 ) 
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Figure 3. Evolution of the spiral arms and the m = 2 mode every 20 Myr. The first column shows the distribution of the stars in polar 
coordinates (i.e. R-4> plane). The peak phases of the m = 2 mode are indicated by dotted lines whose colours (rainbow) correspond 
to the m = 2 mode amplitudes. Time evolution is shown in a rotating frame of the bar with approximately 47 km s“^ kpc“^. The 
horizontal axis indicates the azimuthal angle (p [degrees], and the vertical axis denotes the galactocentric distance R [kpc]. The remaining 
three columns (left to right) plot the amplitude \Bm= 2 \-, pitch angle im =2 [degrees], and pattern speed f^phase,m =2 [km s~^ kpc~^], 
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Figure 4. Evolution of the spiral parameters at each radius (width of 1 kpc) from 2.5 Gyr. The upper horizontal axes shows the time 
normalized by the rotational period Trot at each radius. The spikes in evolution of pitch angle are caused by ill-determined Fourier modes 
when the amplitudes are very small. 
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Figure 5. Bottom: Evolution of the m = 2 modes at /? = 4, 6, 8, and 10 kpc on the im =2 — \Bm =2 \ plane. Each dot represents a value 
sampled every 1 Myr. Top: Amplitude-weighted frequencies of the appearance of pitch angle im =2 at each radius. The hatched region 
corresponds to the predicted maximum pitch angle 2swa around the analysed region due to the swing amplification (see Section iMt . 


where T = —d\nQ/d\n R is a dimensionless shear rate. This 
equation predicts that iswA is a de creasing function of T. 
This trend matches observations (e.g. Seigar_et_aJj 200^ and 


V-body simulations of disc galaxies I Grand et al. 2013li . 


In bottom panels in Fig. [5] we plot the expected values 
of iswA as shaded regions. In the outer regions at i? = 8 
kpc and 10 kpc, the peak pitch angles (around 25°-30°) of 
the simulations match the predictions well. Similar results 
have already been obtain ed for V-body si mulations of mul¬ 
tiple arm spiral galaxies (iBaba et al.ll2013h . In contrast, the 
peak pitch angles at the inner regions (i.e. R = 4 kpc and 
6 kpc) are somewhat larger than the predicted values due 
to the effect of the central bar. Furthermore, we also com¬ 
pare the frequency of appearance for im =2 with iswA- The 
top panels in Fig. [5] show that the most frequent pitch an¬ 
gle matches the values predicted by the swing amplification 
theory, except for R — 4 kpc. 


In addition, the swing amplification theory can explain 


the evolutionary behaviour between r 2 phase,m =2 and \Bm= 2 \ 
because the amplific ation is most ef f icient around the co¬ 
rotating radius (see iDobbs fc Babal |2014| . and references 
therein). In fact, the top panels in Fig. [6] show that the 
most frequent phase speed decreases as a function of the ra¬ 
dius, but it is the same as or slightly faster than the galactic 
angular speed, except for R — 4 kpc. The previous V-body 
simulations of barred spiral galaxies have also pointed out 
that the spiral arms rotate with a decreasing angular rota¬ 
tion s peed that is slightly faster than the galactic angular 
speed llGrand et al.ll2012a^ . However, these studies have not 
been compared with the swing amplification theory. 


3.3 Effects of the bar on spiral dynamics 

In the above subsections, we presented the short-term be¬ 
haviour of the grand-design spiral and showed that these 
results can be explained by the swing amplification mech- 









































































8 J. Baba 


c 

o 

CD 

it 


C'J 


E 

CQ 



Qphase,m=2 [klTl/s/kpc] Qphase,m=2 [km/s/kpc] nphase,m=2 [ktTl/s/kpc] Qphase,m=2 [klTl/s/kpc] 


Figure 6. Bottom: Evolution of the m = 2 modes at = 4, 6, 8, and 10 kpc on |Bm=2|-f2phase,m=2 plane. Each dot represents a 
value sampled every 1 Myr. Top: Amplitude-weighted frequencies of the appearance of phase pattern speed f!phase,m =2 at each radius. 
The hatched region indicates the angular speed of the galaxy Q{R) at each radius. Note that the bar’s pattern speed (at iJ < 3 kpc) is 
approximately 47 km s“^ kpc~^. 


anism. However, this mechanism cannot explain the be¬ 
haviour in the inner regions (e.g. i? < 4 kpc). In this subsec¬ 
tion, to clarify behaviour of the spiral at the inner regions, 
we investigate how the bar affects the dynamics of the grand- 
design spirals. 

To clarify the cone of influence of the bar, we plot the 
results for the (spiral) amplitude \Bm= 2 \ with respect to 
the bar amplitude |i3bar| in Fig. [T] Here, we measure |Hbar| 
with the amplitude of the m = 2 mode at i? = 3 kpc. At 
i? = 4 kpc, there is a clear trend showing that, as the spi¬ 
ral amplitudes become larger, the bar becomes stronger. We 
see a similar trend at iZ = 6, but it has considerable scatter 
around the correlation. In contrast, the results at i? = 8 and 
10 kpc show no clear correlation. These results suggest that 
the bar affects the dynamics of the spiral arm well beyond 
its co-rotation radius (~ 4 kpc; see Fig. and possibly < 6 
kpc. This maximum radius corresponds to approximately 
1.5-2 bar radius and is co nsistent with ob servational val¬ 
ues (approximately 1.4-1.6: ISalo et ^l2010t) . We, therefore, 
suggest that the bar drives the spiral arm within approxi¬ 
mately 1.5-2 bar radius. 

To see the dynamical behaviour of the bar-driven spi¬ 
rals, we investigate whether the inner spiral arms are consis¬ 
tent with the predictions from the invariant manifold theory. 
As described in Section [T] the manifold theory provides two 
predictions: (1) the spirals in strongly barred galaxies will be 
more open than those in less strongly barred ones, and (2) 
the spiral arms should be a bundle of orbits guided by the 
manifolds so that stars should move along the arms rather 
than across them. We compare the simulated barred spirals 
with these predictions. Fig. [8]shows the results for the spiral 
pitch angle im =2 versus |Hbar|- We cannot see any correlation 
between im =2 and |Hbar| for R = 8 and 10 kpc. In contrast, 
for R = 4 and 6 kpc, we see that the low pitch angle spirals 
(*m =2 fy 20°) are lacking for stronger bars (|Hbar| ^ 0.2). 
In this sense, our results correspond to the predictions from 
the invariant manifold theory, although there is considerable 
scatter for weaker bars. Furthermore, to see whether stars 
move along the spiral arm as predicted from the invariant 
manifold theory, we select the star particles located at the 


grand-design spiral at t = 2.78 Gyr and then follow the tra¬ 
jectories of these particles from 2.70 Gyr to 2.84 Gyr. These 
trajectories are shown in Fig. [9] At t < 2.76 Gyr, the chosen 
stars are located at both sides of the arm, and these stars 
form the arm at t ~ 2.78 Gyr. After that, these stars escape 
from the arm to both sides. Therefore, the motions of the 
stars therefore are not simply consistent with the prediction 
fro m the in v ariant manifold theory. 

it rue 3 ll2015h suggested that co-rotating spirals with 
bars are supported by resonant orbits excited in outer Lind- 
blad zones, which are the parents of a continuous sequence of 
resonant radii. To investigate whether this is true, we plot¬ 
ted sample orbits of these particles in the rotating frame in 
Figs. IlOh and llOb . We could not find any resonant closed 
orbits, rather most orbits show monotonic centre shifts of 
the radial oscillations (i.e. radial migrations of stars; Fig. 


galactic discs 

Sellwood & Binneyl 

2002 

iGrand et al.l 

2012a: Baba et al.l 

2013 


Grand et al. 12014 ). This result suggests that the collective 


motion of stars might be an essential ingredient for generat¬ 
ing dynamic spirals due to a non-linear combination among 
the shearing flow, non-linear epicyclic (or p-ellipse) motions, 
and self-gravity (see Section mi). However, this orbit anal¬ 
ysis is limited to the stars selected from a single snapshot. 
Therefore, it would be interesting to further explore the or¬ 
bits of stars comprising dynamic spirals in numerical sim¬ 
ulations for comparison with kinematic theories of stars to 
aid b ar and spiral formation (E omero- G6mej[2012l : IStruckI 
l2015l) . 


4 CONCLUSIONS AND IMPLICATIONS 

In this study, we performed a 3DA-body/SPH simulation of 
a barred spiral galaxy with our origin al A-body/SPH code 
ASURA-2 dSaitoh fc Maki]i3l2009l . [2010l ) and investigated the 
short-term evolution of grand-design spirals in barred galax¬ 
ies. We summarize our main results and give implications 
below. 
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R = 4 kpc R = 6 kpc R = 8 kpc R =10 kpc 






Figure 7. Bar amplitude |-Bbar| with respect to the spiral amplitude \Bm= 2 \ at i? = 4, 6, 8, and 10 kpc. The bar amplitude |-Bbar| is 
measured with the amplitude of the m = 2 mode at R = 3 kpc. Each dot represents a value sampled every 1 Myr. 
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Figure 8. Bar amplitude |-Bbar| spiral pitch angle im =2 at /? = 4, 6, 8, and 10 kpc. Each dot represents a value sampled every 1 
Myr. 
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Figure 9. Motion of selected star particles with respect to the rotating frame of the stellar bar. The star particles are selected from the 
particles associated with a grand-design spiral at t = 2.78 Gyr. Time evolution is shown in a rotating frame of the bar with approximately 
47 km s~^ kpc“^. 
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Figure 10. (a) Sample of the orbits of the stars shown in Fig.[3 Time evolution is shown in a rotating frame of the bar at approximately 
47 km s“^ kpc“^. (b) Time evolution of the radii of the stars shown in panel (a), (c) Probability distribution of the radial migration AR 
[kpc] of the stars shown in Fig. |9] Here, AR = R 2 — Ri, where .Ri and R 2 are the time-averaged radii of stellar orbits at 2.5 < t < 2.8 
Gyr and 3.0 < t < 3.3 Gyr, respectively. 


(i) As suggested b y previous simulations (e.g. 
ISellwood &i Spar^ 1 19881 : iRautiainen fc Sai3 Il999l ~). clear 
grand-design spirals in barred galaxies appear only when 
the spirals connect with the ends of the bar. This implies 
that barred galaxies with spirals not starting from bar ends 
will have weaker spirals than those in barred galaxies with 
spirals starting from the bar ends. Actually, some barred 
spiral galaxies exhibit a clear bar-spiral phase differences 
(e.g. NGC 3124 and NGC 3450: iButa et al.]l2007^ and their 
spirals are not clear grand-design spirals. Yet, this is not a 
statistical argument. To confirm our implication, it requires 
a statistical comparison between barred spirals with and 
without the bar-spiral phase difference. 

(ii) Grand-design spiral arms in a barred galaxy are not 
rigidly rotating patterns but rather differentially rotating 
dynamic pattern s , as suggested by previous simulations 
llBaba et al.ll2009l : lGrand et al.ll2012al l. This means that the 
amplitudes, pitch angles, and phase speeds are not constant, 
but change within a time span of 1-2 rotational periods at 
each radius. Such behaviour suggests that the assumption of 
rigidly rotating grand-design spirals is not necessarily rea¬ 
sonable for investigating gas flows and cloud formation in 
barred spiral galaxies. This point will be disc ussed further 
in Baba, Morokuma-Matsui & Saitoh (|2015bl ') and a forth¬ 
coming paper (Baba et al. in preparation). 

(iii) The stellar bar affects the evolution of spiral arms 
within 1.5-2 bar radius, and this ra dius is consistent with ob¬ 
servations of barred spiral galaxies (|Salo et al.ll201fll ~l. Within 
this radius, as predicted from t he invariant manifold theory 
fe.g. lAthanassoula et al.ll2009bl ~l. low pitch angle spirals are 
lacking for stronger bars, although there is considerable scat¬ 
ter for weaker bars in the bar amplitude-spiral pitch angle 
plane. In addition, motions of stars located at a grand-design 
spiral are not simply consistent with the prediction from the 
invariant manifold theory. 

(iv) The swing amplification theory can explain the most 
frequent pitch angles and the amplitude-pitch angle rela¬ 
tionships in the outer regions (i.e. R > 1.5-2 bar radius). 
These results imply that interplay among the galactic shear¬ 
ing motion, epicyclic motion of stars, and self-gravity of stars 
is an essential factor for the dynamics of spiral arms even in 


barred spiral galaxies. Note that we assumed that the DM 
halo is a fixed external potential, so th at the bar is formed 
via bar instabilities in ma ssive discs (lOstriker fc PeeblesI 
I 1973 I : iMstathiou et al.lll982h . In this case, the self-gravity of 
the disc plays a role in the formation of spiral arms via the 
swing amplihcation. Gonsidering a ‘live’ DM halo, however, 
lAth anassoui^ (l2002l l showed that bars can develop even in 
less massive discs. We infer, then, that the spiral arms in 
such discs may show different behaviour from those in mas¬ 
sive discs such as those investigated in this study. This topic 
will be the subject of future work. 
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